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The  number  of  test  samples  used  to  characterize  the  fatigue  constants  needed  for 
failure  predictions  for  ceramic  materials  determines  the  confidence  in  these 
predictions.  The  statistical  reproducibility  of  the  dynamic  and  static  fatigue 
experiments  used  to  measure  the  fatigue  constants  was  analyzed  using  both  sta¬ 
tistical  theory  and  a  fionte  Carlo  computer  simulation  technique.  It  was  found 
that  the  statistical  reproducibility  depended  not  only  on  the  number  of  test 
samples  but  also  on  the  other  experimental  test  variables.  It  v.'as  (continued) 
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Block  20.  Abstract  (continued) 

shown  that  the  uncertainty  in  the  statistical  reoroduci bi 1 i ty  can  be  large 
especially  for  sample  size  less  than  about  IOC.  Guidelines  for  selecting 
the  optimun'  samole  size  for  a  giver,  ayrarric  O’"  static  fatigue  experiment  are 
given.  It  is  recommended  that  before  meaningful  conclusions  can  be  drav.'n 
regarding  the  effect  of  a  test  variable  on  ^atigue,  the  statistical  rep’"oauc- 
ipility  of  tne  experiment  be  determined. 


FOREWORD 


This  report  describes  the  results  of  an  experimental  program 
oriented  tov/ard  a  better  understanding  of  lifetime  predictions  for 
optical  glass  fibers.  Some  of  the  progress  made  toward  this  goal 
is  summarized  in  the  attached  technical  paper  comprising  this 
report. 
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ABSTRACT 

The  number  of  test  samples  used  to  characterize  the  fatigue  constcrt: 
needed  for  failure  predictions  for  ceramic  materials  deterndnes  the  cori- 
fidence  in  these  predictions.  The  statistical  reproducibility  of  tne 
dynamic  and  static  fatigue  experiments  used  to  measure  the  fatigue 
constants  was  analyzed  using  both  statistical  theory  and  a  Monte  Carlo 
computer  simulation  technique.  It  was  found  that  the  statistical  reproduci¬ 
bility  depended  not  only  on  the  number  of  test  samples  but  also  on  the  other 
experimental  test  variables.  It  was  shown  that  the  uncerta i nt ■/  in  the  statistical 
reproducibility  can  be  large  especially  for  sample  size  less  than  about 
100.  Guidelines  for  selecting  the  optimum  sample  size  for  a  given  dynamic 
or  static  fatigue  experiment  are  given.  It  is  recommended  that  before 
meaningful  conclusions  can  be  drawn  regarding  the  effect  of  a  test 
variable  on  fatigue,  the  statistical  reproducibility  of  the  experiment 


be  determined. 


1.0  INT;'’0DU:TI0N 


Methods  dealing  with  design  protlems  inv'ol  vine  fetioue  of  Cfcran.ic 
materieiS  neve  been  developed  over  the  past  10  years  throogr.  tne  applicatior' 
of  fracture  iiiecnanics  principles.  Since  these  principles  car  oe  used  to 
characterite  both  Lhe  conditions  for  suDcritical  crack  growth  and  the 
conditions  for  crack  instability,  they  can  be  used  for  purposes  of  design 
to  estimate  the  allowable  stress^or  the  expected  lifetime,  or  tne  proc^ 
stress  necessary  to  assure  a  minimum  lifetime.  This  is  accoir.pl  i  sriec  by 
estimating  the  initial  crack  size  in  a  ceramic  component  end  tne  time 
required  for  this  initial  crack  to  grow  to  a  critical  size  for  sbontaneous 
fractu’“e.  For  example,  it  has  been  derived  by  assuming  a  simple  ,nower  law 
relationship  between  subcritical  crack  velocity  and  stress  intensity  tnat 
the  failure  time  (t,)  under  a  constant  applied  stress  (c  )  is;^’'"^ 

T  S 


■)-  r  B 


/  2  N-2 

where  B  =  2/(AY  (N-2)  K.^  ),  A,  N  =  material /environment  constants,  Y  = 

geometric  constant  (about  1.2  for  surface  flav/s),  =  critical  stress 
intensity  facto”,  and  '  fracture  strength  in  an  inert  environment.  B  and 
N  in  Eq.  (1)  are  fatigue  constants  that  for  a  given  material /environment  system 
characterize  subcritical  crack  growtn.  The  inert  strength  in  Eq.  (1)  character¬ 
izes  the  initial  flaw  size.  If  p'oof  testina  is  used  to  truncate  the  flav; 
distribution,  then  the  miinimum  inert  strength  after  proof  testing  is  equal 
to  the  maximum  proof  stress  ("t^),  hence,  tne  correspondingly  minimum  failure  time 


^  -  a. 


Froiii  Lq.  ■]  citicJ  (2';  Tt  ^  s  seen  that  failure  i.red  u  nntii  art-  nMit-riOt-r  •, 
ori  tiie  Tatiquij  iiaraiiieter:.  N  and  ii.  Trf’Se  t'J '  art;  v  tm-.  laii'.  ■.  'Di  c. 

given  ladteria  !  ,.'env  1  roriiiient  system  and  can  be  experimentally  duti-r;:  uu-d 

3 

directly  using  fracture  mechanics  techniques,  or  can  be  inCirectl>  measured 

2 

using  static  fatigue  or  dynannc  fatigue  experiments.  unfortunate','.-,  foiiure 
predictions  are  extremely  sensitive  to  the  experimental  unc  erta  i  nt m  trie 
fatigue  parameters.  Statistical  techniques  for  estimating  trns  uncertair,t\ 
in  failure  predictions  have  been  developed;  '  however,  the  statistica’ 
reproducibility  of  the  experimental  techniques  used  to  evalu.ne  itie  fatigue 
parameters  !J  and  5  has  not  been  previously  determined.  Statistical  repro¬ 
ducibility  is  due  to  random  sampling  errors  that  are  inherent  in  every  ex¬ 
periment.  In  the  random  selection  of  a  finite  number  of  samples  for  testing, 
one  would  expuct  to  see  some  statistical  variability  in  the  measured  properties 
that  would  be  dependent  on  the  number  of  samples  selected.  It  further  would 
be  expected  that  this  variability  in  the  estimation  of  the  fatinue  oaraineters 
would  increase  as  sample  size  decreases;  however,  this  statistical  reproduc i bi 1 i 
has  never  been  quantified  although  some  Monte  Carlo  results  on  the  re¬ 
producibility  of  N  as  determined  by  dynamic  fatigue  tests  have  been  previously 
reported  by  the  present  authors.^  Before  meaningful  conclusions  can  be  drawr- 
from  the  results  of  fatigue  experiments,  the  statistical  reproduc i b i 1 i ty  of 
these  experiments  must  be  known. 

The  purpose  of  this  paper  is  to  quantify  the  statistical  reproducibility 
of  the  dynamic  and  static  fatigue  experiments  for  measuring  the  fatigue 
constants  K  and  E.  Statistical  reproducibility  is  analyzed  usino  both 


statistical  theory  and  a  Monte  Carlo  computer  simulation  technique.  Since 


the  statistical  theory  contains  a  number  of  critical  assumotions,  it  is 
important  zc  i noepenoently  validate  the  statistical  theory  aDoroach  with 
the  Monte  Carlo  technique.  Emphasis  is  placed  on  tne  dynamic  and  static 
fatigue  tecnniques  because  they  are  increasingly  being  used  to  measure 
the  fatigue  behavior  of  ceramics.  This  is  because  these  test  techniques 
can  utilise  samples  containing  flaws  representative  of  those  on  actual 
components.  It  is  believed  that  the  results  of  this  study  will  leac  to  a 
better  understanding  of  the  static  and  dynamic  fatigue  test  technioues  and 
their  statistical  variability.  From  this  information  guidelines  for  sample 
size  reauirements  in  terms  of  optimum  statistical  reproducibility  can  be 
developed. 


I 


2.0  ANALYSIS  OF  STATISTICAL  REPRODUCIEIirV 


2.1  Dynamic  Fatigue 

Dynamic  fatigue  data  is  generated  by  measuring  the  fracture  strength  o* 
a  number  of  samoles  at  several  constant  stressing  rates.  The  fatigue  constants 
N  and  E  can  be  determined  from  dynamic  fatigue  data  througn  using  one  of  four 
analyses:  median,  homologous  stress,  iterative  bivariant,  or  iterative  tri- 
variant. Since  all  of  these  techniques  analyze  the  same  set  of  dynamic 
fatigue  data  in  determining  N  and  E  through  a  linear  regression  analysis,  all 
are  expected  to  result  in  essentially  the  same  statistical  reproducibility 
for  N  and  B.  Thus,  the  median  analysis  technique  was  chosen  for  this  study 
because  of  its  simplicity  and  wide  usage. 

2  7 

With  the  median  analysis  the  dynamic  fatigue  data  are  fitted  to:  ’ 


f. 


(3) 


where  S  =  median  fatigue  fracture  strength,  c  =  stressing  rate,  and  a, ,  a^,  - 

a  2 

linear  regression  constants.  The  fatigue  constants  are  then  determined  from:  ’ 


>^,  -  N  i  i  )  a. .  -  +  \  _  f  fj,  y  ',  , 

"  "f  I  ■"  \  ~r  ■  3  ;  /'vrvS- 

where  =  median  inert  strength. 

Dynamic  fatigue  strength  data  of  an  "ideal"  material  was  simulated 
on  a  computer  using  a  Monte  Carlo  technique.  It  was  assumed  that  the 


(4a) 

(4b) 


fatigue  constants  N  and  B  of  this  ideal  material  are  given  and  that  the 
inert  strength  distribution  is  given  by  a  two  parameter  Weibull  distribution 


whose  slope  and  scale  parameters,  m^and  respectively,  are  known.  IJith 
the  Monte  Carlo  technique  a  given  number  of  samples  at  a  specT^ic  stress¬ 
ing  rate  were  chosen  by  randomly  selecting  their  fBilu*-e  C'C-bai; i ' " t\ 
from  a  uniforii  distribution  between  C  and  1.  The  correspond i nc  *'atigue 

fracture  strengchs  were  then  calculated  according  to  the  fracture  mechanics 
2 

relationship: 

r  C  --  -  ^  ^  ,'^z  ^  -o.)  *  r  1 

A/./"  • 

A  separate  set  of  inert  strength  samples  were  chosen  similai'ly  to  the  fatigue 

samples,  namely  failure  probabilities  were  ranoomly  chosen  and  the  corresooncing 

1  0 

inert  strength  were  calculated  from,  the  two  parameter  Weibuli  distribution: 


Once  a  set  of  inert  strength  and  fatigue  strengths  at  several  different 
stressing  rates  were  randomly  chosen,  N  and  B  were  determined  from  this  data 
using  Eqs.  (3)  and  (4).  It  should  be  noted  that  the  median  strength  for  a 

given  set  of  strength  values  was  determined  by  ranking  and  fitting  the  strengths 
to  a  Weibull  distribution  by  linear  least  square  method.  The  meaiar  strength 

was  then  calculated  as  the  value  at  F=  0.50.  Alternatively,  tne  median 
strength  could  have  been  determined  by  choosing  the  actual  median  strength 
value;  however,  this  results  in  greater  variability  and  thus  was  not  used.  By 
iterating  this  procedure  100  times,  distributions  for  N  and  B  were  generated 
which  represent  the  statistical  reproducibility  of  N  and  B  as  determined 
from  the  dynamic  fatigue  test.  With  the  100  values  generated  by  the  Monte 
Carlo  technique,  the  average  values  with  their  corresponding  variances  and  the 
covariance  between  N  and  B  were  calculated  from  the  usual  statistical  formulas:''"' 


/J  ' 


^  O  "  n 


\J  f/^)= 


l/(A&]  -  -  I  -  ->’-'3  ) 


C.(wl6i  ^  -  I  i 

^  -nw  j,, 


v/here  n  =  100  for  these  computer  calculations.  For  a  given  materi al /envi ronment  system, 

i.e.  for  a  given  set  of  m,  S  ,  r>i,  and  B,  the  important  test  variables  studied  were 

0 

the  number  of  samples  per  stressing  rate,  the  stressing  rate  range,  and  the 
number  of  stressing  rates  that  were  used  in  a  given  test.  For  the  case  of 
mul ti -stressing  rates,  the  stressing  rates  were  evenly  spaced  from  maximum 
to  minimum.  Figure  1  gives  a  schematic  flow  diagram  of  the  Monte  Carlo 
computer  simulation  technique  for  determining  the  fatigue  constants  N  and  B 
by  the  dynamic  fatigue  test.  Note  that  the  same  number  of  inert  samples  were 
chosen  as  the  number  of  samples  per  stressing  rate. 


1  he  variances  and  covariance  of  the  inert  strencf:  oerar.eters  we“'e  a' 
oerivec  as  a  function  or'  samole  size  fror  statistical  tneo>"v  rrafirr 


certain  assun'.ttions  (see  Apoendix,  Eqs.  (A15(,  (A2c0,  anc  ( 


He  j  <  .  ITic  ^651. 


of  this  analysis  are  as  'ollows; 
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I . "  /  ■ 


;  :  *  1-^  .  ) 


where  J  =  number  of  samples  and  f  =  gai.in'.a  r'unction.  Lii;ewise,  the  var'ances 


the  fatigue  constants  N  and  E  and  their  covariance  were  derivec  as 


s  nown  i n  t 


ppendix,  Eqs.  (A41),  (A42),  and  (A431,  to  be: 


V  'jy.  - 


}<  ,  ;  /J  -  1;  ' 


,  0 .f.  I  -  •  \  -  t‘ 


Cc.jnj  t  - 
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ir  E  ( /V-# \  -  2  uV,  5.  -2  A-  X.  c 
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r. . r  1 


where  =  total  number  of  samples  used  in_  dynamic  fatigue  experimeni 
number  of  stressing  rates,  RiS^  r  z  ~  X  ( £•,  C.  -  i-n  q-  )  /  4 


r.  '  \  iJ  -  ^  j 
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Since  these  de'^ivations  rriade  e  nuiriDer  cf  critical  assumptions,  it  is  o"' 
importance  to  compare  these  variances  witn  tnose  aeneretec  by  tne  lionie  Ca*'lc 
tecnnique  tc  oeterrrine  me  validity  of  che  assumptions. 

2.2  Scatic  Fatigue 

Static  -'atigue  tests  entail  the  repeatec  measurement  o*  *a''".u''e  time  at 

several  constant  appliec  stresses.  Static  fatigue  data  can  oe  analyzed 

similar  to  dynamic  fatigue  data  by  the  median,  homoloaous  stress,  iterative 

2  7 

bivariant,  and  iterative  trivariant  analyses.  ’  As  in  the  case  cynamic 
fatigue  the  median  analysis  was  cnoser,  for  stuoying  the  statistical  repro- 
ducibilitv.  w'itn  tne  meoian  analvsis  the  data  is  fitted  tc; 


t  ^ 


(10) 


where  t.  =  median  failure  time  and  a.,  a,  =  linear  recression  constants, 
f  j  4  ' 

2 

The  fatigue  constants  N  and  B  are  determined  from:  ’ 


(lie) 


(11b) 


The  Monte  Carlo  analysis  of  the  statistical  reproducibility  of  A  and  B 

as  determined  from  static  fatigue  data  is  similar  to  that  used  for  dynamic 

fatigue  data.  First,  a  given  number  of  failure  times  at  a  specific  applied 

stress  are  selected  by  randomly  choosing  the  failure  probability  and  then 

2 

calculating  the  corresponding  failure  time  from; 


A' 


r 


(1 


CsJ 


Ofice  a  set  of  failure  times  at  several  di'^ferent  appliec  stresses  ario  e  set 
cf  inert  strengths  are  ranooir.ly  chosen,  N  and  E  are  determinea  fro;:,  Ecs.  ’10' 


1 


and  (11;.  By  iterafinc  this  procedure  IOC  tiines,  a  distribution  'n  and  E 
values  ere  generated  frorri  which  the  statistical  reoroduci  bi  1 ty  o*  these 
paraiTieters  can  be  aeterir.i ned .  Similar  to  the  case  of  dynamic  *aticjje,  the 
variables  studied  for  a  given  material/environment  system  were  the  number 
samples  per  anplied  stress,  the  aonlied  stress  range,  and  the  nui-aer 
applied  stresses  in  a  given  test. 

The  variances  in  the  fatigue  constants  N  and  E  and  their  covariance  co,.lc 
also  be  derived  from  statistical  theory  as  shown  in  the  Aooencix  ,  Eds.  (A5E  ) , 


(A53),  and  (A54)  to  be; 


I,/;/). 
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(13b) 
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where  now  J  =  number  of  samples  used  in  static  fatiaue  experiment,  =  number 
0  -  a 

of  applied  stresses.  J.|  =  number  of  samples  tested  at  each_applied  stress. 
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tr.,.  >  2,  -  1.^^. 

Again  it  wes  of  inte'^est  to  conipere  these  variances  to  those  dete*'rri ned  by  tne 
Monte  Carlo  technique  to  determine  the  validity  of  the  assumptions  made  in  ceri' 
ing  the  above  equations. 


3.C  RESULTS  AKD  DISCUSSION 


Both  the  Weibull  parameters  (m  and  S^)  and  the  fatipje  constants  (!<  enc 
depend  stroncly  on  samcle  size.  Figures  2  and  3  show  ihe  oependency  *or  n  an 
S  .  resoecci  vel  V  on  semole  size  where  the  coefficient  cf"  variation  is  de'f'inec 

O  ' 

as  the  standard  deviation  divided  by  the  mean  value  of  tne  oarameter  in 
question.  From  tnese  figures  it  can  be  seen  that  tnere  is  good  agreement 
between  the  Monte  Carlo  technique  for  estimating  the  statistical  reprodu¬ 
cibility  and  that  derived  from  statistical  theory.  The  figures  further 
show  that  the  statistical  variability  in  the  Weibull  slope  parameter  m  is 
a  function  of  only  sample  size  while,  for  the  Weibull  scale  parameter 
the  variability  is  a  function  of  both  sample  size  and  the  Weibull  slope. 

Since  for  small  sample  sizes  the  statistical  uncertainty  in  the  Weibull 
parameters  can  be  large,  especially  for  m,  important  judgements  and 
significant  analyses  of  strength  should  not  be  based  on  small  sample  sizes. 
From  Fig.  2  ana  3  it  can  be  seen  that  sample  sizes  of  at  least  30  should  be 
used  for  all  but  the  most  preliminary  investigations,  although  for  small  m's 


to  get  acceptable  levels  of  5^  may  require  as  many  as  100  samples.  Similar 
results  and  conclusions  were  reached  in  earlier  Monte  Carlo  simulation 


studies . 


12,  13 


Figures  4  and  5  show  that  the  statistical  reproducibility  of  the  fatigue 
parameter  N  and  B  as  determined  by  the  dynamic  fatigue  experiment  is  strongly 
dependent  on  both  sample  size  and  the  fatigue  resistance  of  the  materia' 
(large  N  values  generally  represent  materials  with  a  greater  fatigue 
resistance).  For  sample  sizes  of  less  than  100,  the  statistical  uncertainty 


in  N  and  B  can  be  very  large,  especially  for  the  more  fatigue  resistant 
material.  Again  there  is  good  agreement  between  the  Monte  Carlo  technioue  fo 


theory;  thus,  giving  evidence  of  tne  validity  of  the  assumotions  rr.dde  in 
deriving  the  equations  for  statistical  variability  (see  Appendix). 

Tne  statistical  variability  of  the  fatigue  parameters  !i  arid  B  as  determine 
by  tne  oynamic  fatigue  experiment  was  also  dependent  on  the  nurribe*"  ci"  stressing 
rates  chosen,  stressing  rate  range,  and  the  Weibull  slope  rr,.  For  examrle. 

Fig.  6  shows  that  for  the  same  range  of  stressing  rates  (maximum  to  minimum) 
and  the  same  total  number  of  samples,  uncertainty  in  statistical  reproduci bi 1 i t 
increased  as  the  number  of  stressing  rates  used  in  determining  K  is  increased 
from  2  to  7.  The  best  statistical  reproducibil ity  occurs  for  the  case  where 
N  is  determined  from  strength  measurements  at  two  stressing  rates  corresoondi nc 
to  the  maximum  and  minimum.  Figure  7  shows  the  statistical  variability  of 
as  a  function  of  m,  keeping  the  other  oarameters  constant,  and  illustrates 
that  low  m  values,  corresponding  to  a  greater  variability  in  strength,  re¬ 
sult  in  a  larger  uncertainty  in  fl.  Figure  8  shows  that  the  statistical 
variability  of  II  is  quite  sensitive  to  the  stressing  rate  range,  with  the 
variability  increasing  with  a  decreasing  range  of  stressing  rates.  Finally, 
it  should  be  noted  that  the  Monte  Carlo  results  were  left  out  of  Figs.  6,  7, 
and  8  for  clarity;  however,  these  results  agreed  quite  well  with  those  shown  in 
the  figures  based  on  statistical  theory. 

The  distributions  of  the  fatigue  parameters  N  and  B  as  generated  by  the 

Monte  Carlo  computer  simulation  technique  for  the  dynamic  fatigue  experiment 

could  be  approximated  by  a  normal  distribution.  Finure  ?  shows  that  a  normal 

distribution  well  represents  the  histogram  for  the  fatigue  parameter  h. 

Histoarams  for  the  Weibull  parameters  m  and  S  could  also  be  approximated  bv 

0 

normal  distributions.  This  is  important  because  the  confidence  limits  of 


reproducibility  for  a  given  parameter  can  be  estimated  by  simply  multiplying 


ne  stancara  aev’aaion  bv  bne  accrocr': aaa  ‘actor  -cr  nc'-'ia.  :i3:r‘c;. 


-115,  cna  reorcGcci  ci  1 1  ty  l-nits  wcu'c  be  .  T-cC  acapcarc  cev  -  a  cvcrs  '  arc 


'-ecresent  cce  'ncarvai  ,‘<nera  aocut  rbb  tr  ere  va.'jes  "cr  a  :i /en  aa.'-aceca''  .’e. 
~be  rCb  reorcccc:  bTi '  cy  ’'ciica  .vculc  ce  star, care  cevac-'srs  '  arc  sc  an. 

Rescics  :r  -re  statisccca'  recrcCLici b- ' : ty  C'  ere  'ac^cce  car3;~e:ars 
.'i  ana  3  as  cetarminea  by  the  static  facicue  exceri'ent  were  sirci "!  ar  cc 
-.resa  s.rcwn  -or  tre  dynarnic  fatigue  experimenc.  In  paraicuiar,  tre  statistical 
7arsa3i’:ty  of  d  and  3  as  beteninea  oy  the  static  -‘atigue  excerc”er,t  cecreasec 
with  increasing  samole  size,  deC'^easad  for  nater-als  'w-itr  ^ovier  'i  vaiues, 
tecreasec  as  tre  acoliec  stress  '■ance  increased,  and  decreasec  as  the 
weipuii  sloce  ca  •ar.ate’-  r  ■'ncreasec.  -or  the  same  aorrec  stress  range  arc 
numce'-  b*  total  samoles,  tr.e  best  reorccuci bi  1  •  ty  occurs  for  the  case  wrere 
.!  ana  3  are  cetermi nee  -rem  time-tb-’^a’ ■  ure  'neasurements  at  two  aboI''ec 


isses  corresbonci  rg  to  tre  'naximuin  anc  minimLirr,.  For  one  case  :-  ~ul  t' -aoo' •' ec 


stresses,  the  statistical  variability  depended  nainiy  on  the  total  number  of 
samples,  not  on  the  numoer  of  samples  per  applied  stress,  and  only  somewnat 
on  the  number  of  applied  stresses  used  to  determine  N  and  3.  As  before. 


there  was  gocc  agreement  petween  the  .Venta  Carlo  analysis  and  the  statistically 


cerivea  vanaoiiity  ecuations  .see  .-ippenc-x ; .  Aiso,  the  ci  stnoutions  generatec 


by  the  Nionta  Carlo  analysis  could  be  approximated  by  normal  distributions. 
Finally,  it  should  be  noted  tnat  the  magnitude  of  the  statistical  reprocucib- 
of  ‘i  and  5  dete.'-mihed  by  tna  stat'c  -'atigue  excer-'ment  is  svmi'ar-  to  tnat 
oeterm'hec  -rpm  the  evra.m-o  -aticue  exoer-ment  -or  tveioa'  '■anoas  of  o'/nj-',: 


.V 


1 


h 
\  ■ 


i 
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Finally,  it  should  be  noted  that  there  was  a  small  bias  present  in  the 
Monte  Carlo  generated  values  for  ni,  S^,  N,  and  B,  especially  for  small 
sample  sizes.  It  is  well  known^ that  when  random  variables  are  combined  in  a 
non-linear  fashion,  the  resulting  quantity  is  generally  subject  to  biases,  i.e. 
the  combined  effect  of  the  random  fluctuations  of  the  individual  variables 
will  cause  the  derived  quantity  to  be  systematically  larger  or  smaller  than 
it  would  have  been  in  the  total  absence  of  such  fluctuations.  Biases  in  m 

1  O  10  1  ^ 

and  have  been  previously  discussed.  ’  ^  These  systematic  biases  in  the 

derived  parameters  were  quite  small  in  comparison  with  the  statistical 
reproducibility  of  the  parameters;  hence,  they  were  neglected  in  the 
statistical  reproducibility  of  a  given  parameter. 


.  0  . 


^r;c~:cal  :::!s::£RAT:o:is 

Once  “ne  sta^'is'ica'  ''eoracuc' bi  ’  i  ty  3~  tne  'aaiaue  :ararr,aae''3  'i  ara  £ 

IS  '<r.ovvn  as  a  'uncaian  S'  sa~c’,e  s'za  anc  ■r.a  atner  aess  /ar'ao'es,  a.-'e 

carresacr.ci nc  saaa’sa'.aa'  jaccraa:".ay  in  -ai’.jre  araciaaians  :an  aaer  :e 

1  “ 

aeaar-'.i  nec  as  aisajssaa  :''a7':us'y.”’'  ::  is  -^Ofa  anis  knew' acre  anaa  ane 

nuhicer  of  samo'es  racL.’re'0  *or  a  given  oynamic  or  saaaic  'aaigue  aasa  :an 
be  deaenminec.  ~'ne  e.xaoa  s.'^c'ce  a*  san'c'a  size  *ar  a  aiven  aesa  wi’’  oeoe'^d 
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in  the  dynamic  fatigue  test.  An  uncertainty  of  +  lO'i;  requires  about  140 
samples  to  be  used  in  the  dynamic  fatigue  test  for  the  stressing  rate  range  of 
0.013  to  50  MPa/s.  However,  it  can  be  seen  from  the  Fig.  12  that  if  the  stress¬ 
ing  rate  range  is  increased  to  0.00G5  to  50  .MPa/s,  then  the  number  of  samples 
required  woulo  decrease  to  about  75  On  the  other  hand,  if  the  stressing 
rate  range  is  decreased  to  0.05  to  5.0  MPa/s,  the  number  or  s.amples  required 
increases  to  about  300.  This  figure  clearly  shows  that  sample  size  require¬ 
ments  are  not  necessarily  small  and  depend  on  the  specific  degree  of  reproduci¬ 
bility  that  is  acceptable. 

The  analysis  of  statistical  reoroduci bi 1 i ty  can  also  be  useful  in  de¬ 
termining  whether  a  given  variable  effects  the  fatigue  benavior  of  a 
material.  For  example,  the  fatigue  behavior  of  soda-lime  glass  has  been 
measured  as  a  function  of  test  environment  (6N  NaOH,  distilled  water,  and 
6N  HCl)  using  the  dynamic  fatigue  test  with  stressing  rate  range  of  .;6  to 
0.17  MPa/s  and  a  total  of  120  samples  for  each  dynamic  fatigue  test.^'^ 

The  effect  of  test  environment  can  be  seen  by  calculating  the  allowable  stress 
for  samples  that  have  been  proof  tested  up  to  ICO  MPa  and  that  must  survive 
a  minimum  of  one  year  in  service.  The  result  of  these  calculations  are  given 
in  Table  I  with  the  indicated  +  one  standard  deviation  limits  due  to  the 
statistical  uncertainty  in  the  fatigue  constants  N  and  B.  It  is  seen  that 
the  three  predicted  allowable  stresses  are  well  outside  of  a  standard  deviation 
from  each  other  and,  thus,  it  is  likely  that  the  test  environments  do  have  a 
significant  effect  on  the  fatigue  behavior  of  soda-lime  glass.  While  this 
conclusion  seems  quite  straightforward,  it  must  be  remembered  that  before  any 
meaningful  conclusions  can  be  drawn  regarding  the  effect  of  a  particular 
variable  on  the  fatigue  behavior  of  a  material,  it  must  be  demonstrated  tnat 
the  effect  is  larger  than  the  statistical  reproducibility  of  the  experiment. 
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Table  1.  Predicted  Allowable  Applied  Stress  for  a  Minimum  Lifetime  of 
One-Year  for  Soda  Lime  Glass  Samples  Proof  Tested  at  100  MPa 


Environment  _ N _ 1  n  3  (MPa^ ■  s  )  "a 


6N 

NaOH 

19. 

,5 

(+  2.81) 

0. 

,163 

it 

1.64) 

25. 

,95 

(l 

3. 

.36) 

Dis 

tilled  H^O 

13, 

.0 

(+  1.11) 

2 , 

,585 

(  + 

0.87) 

1  5 . 

,92 

(t. 

1  . 

,53) 

6N 

HCl 

25. 

,1 

(+  4.35) 

-16. 

,058 

(1 

4.67) 

13, 

.36 

(i 

2 , 

.15) 

+  Number  in  parent.hesis  is  +  one  standard  deviation  corresponding  to  statistical 
reproducibility  of  dynamic  fatigue  experiment. 
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T'ne  coefficienL  of  variation  (C.V.)  of  Weibull  sloce  oaramete'"  m 
as  a  function  of  sample  size  ^or  various  values  of  m.  The  Tcnte 
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predictions  from  statistical  theory,  Eg. (3a;. 

The  coefficient  of  variation  of  Weibull  scale  parameter  S  for 
various  values  of  m  as  a  function  of  samole  size.  The  Moflte  Carlo 
results  are  given  by  the  data  ooints  and  are  compared  to  crec'cticns 
from  statistical  theory,  Ea.(3b). 

The  coefficient  of  variation  of  fatigue  parameter  'i  as  a  function 
of  total  number  of  samples  used  in  dynamic  fatigue  experiment  .vith 
3  stressing  rates  (0.013  to  50  MPa/s).  The  Monte  Carlo  results  are 
given  by  the  data  ooints  and  are  compared  to  predictions  from  statisti¬ 
cal  theory,  Ea.(3a;.  Th^  '.'ieibull  and  fatigue  constants  chosen  are 
typical  of  soda-lime  glass'  (.m-=  15.0,  S  .=  133.0  MPa,  3  =  13.4,  and 

?  '  ^  "I  ■'  J 

B  =  0.13  MPa“‘s)  and  vitrified  grincinc  wheel  material' 

(;7i^=  15,  5^.=  53  MPa,  M  =  <-3.2  and  3  =  1.1  x  10'^  >'Pa'-s’, 

The  coefficient  of  variation  of  fatigue  parameter  3  as  a  ‘unction  of 

total  number  o’"'  samples  used  in  dynamic  fatigue  experiment  witn  3 

stressing  rates  ^0.013  to  50  MPa/s).  The  Monte  Carlo  results  are 

given  oy  tne  cata  points  and  comcared  to  precicticns  ‘rom  statistical 
theory,  Ec.lPbl.  The  Weibull  and  fatigue  constants  chosen  are  typical 

of  soda-lime  giass^  (m.=  13.0,3^.=  133.0  MPa,  N  =  13. d,  and  3  = 

9  ,  1  01  ■,  j 

0.13  MPa'-5i  and  vitrified  grinding  ./heel  material'  IE, 

3  .=  53  '-’Pa,  .'1  =  43.2  and  3  =  1  . 1  x  lO'"^  MFa'-s). 
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The  coefficient  of  variation  of  fatigue  parameter  '!  as  a  ‘unction  o* 

total  number  pf  samples  used  in  aynamic  ’"atigue  experiment  ■■/itn 

stressing  rate  range  of  0.313  to  50  MPa/s.  The  number  of  stressing 
rates  used  in  determining  ■'1  is  as  indicated  and  the  appropriate  constants 

are  11  =  13.4,  3  =  g.13  \ipa~-s,  m.=  3.2,  and  S  .=  133.0  MPa. 

1  0 1 

'he  coe'ficient  of  variation  of  fatigue  carameten  h  as  a  ‘unction  of 
total  number  of  samples  used  in  dynamic  ‘‘aticue  expe'"'i'en‘  3 
stressing  ^ates  (3.0:3  to  50  MPay-s)  aria  as  a  ‘unction  oh  Te-bull 
5 1  ooe  parameter  while  tne  otne'’  cpnst3'’ts  a^'e  'I  =  13.4 

3  =  0.13  MPa~’S,  ana  5^--=  130.0  MPa'. 


!^iGure  3.  The  coefficient  of  variation  of  fatigue  paraireter  T  as  a  'jnctian  of 
total  number  of  samples  useo  in  dynamic  fatigue  exoeriment  'or  3 
stressing  rates  with  two  ranges:0.013  to  50  ana  0.35  to  5.0  ^f'd/S. 

The  appropriate  constants  are  'I  =  13. 'i,  3  =  0.13  MPa'-s,  m.=  5.2, 
ana  Soi=  138.0  M.°a. 

Figure  9.  .histogram  for  fatigue  parameter  11  as  cetermined  Monte  Carlo  si'^ulated 
dynamic  fatigue  experiment  for  3  stressing  rates  yO.C13  to  50  ’-'Pa,  s' 
using  a  total  of  150  samoles.  The  supermooseo  smooth  curve  is  the 
probability  density  function  of  the  normal  distribution  witn  t.-'e  same 
mean  and  variance  as  the  'lonte  Carlo  results. 

Figure  10.  The  coefficient  of  fatigue  parameter  'I  as  a  function  of  total  number 
of  samples  used  in  dynamic  fatigue  exoeriment  with  3  stressing  ■"ates 
(0.013  to  50  MPa/s)  and  static  fatigue  exaeriment  with  3  aoolied  stresses 

(40  to  53.6  MPa).  Experiments  are  tyoical  for  soda-lime  glass" 

(N  =  13.4,  B  =  G.13  MPa^-s,  m .  =  3.2,  and  S^.=  13.S.0  MPa). 

Figure  11.  Design  diagram  for  minimum  lifetime  prediction  after  orcof  testing  i^or 

soda  lime  glass  (M  =  13.4,  3  =  0.18  MPa'-s,  m, =  3.2,  and  S  .=  133.3  ''Pa) 

I  0 1 

The  90t  confidence  limits  are  shown  for  3  samole  sizes  (45,  30,  and  13C) 
used  in  dynamic  fatigue  experiment  with  5  stressing  rates  (0.013  to 
50  MPa/s). 

"^igure  12.  90“i  confidence  limits  for  the  proof  stress  to  assure  a  minimum  lifetime 
of  10  years  at  a  constant  aoplied  stress  of  7  MPa  for  soda-lime  glass 

(11  =  13.4,  3  =  0.18  MPa^-s,  m;=  3.2,  S  .=  138.0  MPa)  as  a  -'unction 
of  sample  size  used  in  dynamic  fatigue  e.9berirr,ent  with  3  stressing  rates 
in  the  range  of  0.05  to  5.0,  0.013  to  5C,  and  0.005  to  50  '1Pa/s. 
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APPENDIX 


A.'.  Linecr  Regression 

Linear  regression  is  gene"al'!y  usee  zc  analyze  tne  test  data  in  oete'-"'’' n- 

inc  one  fatigue  oaranieters  N  and  E.  Tne  variances  o*  these  ■fat'aue  oe’'a''iete’‘S 

can  Chen  be  estimacec  dv  using  cne  lai\  C"^  orooagation  of  e'^rors.  commonly 

11  ' 

known  as  the  chain  ruie.  ,ne  chain  rule  is  a  linear  approximation  and 

becomes  inaccurate  when  the  coefficient  of  variation  [ratio  of  the  standard 
deviation  tc  tne  mean'  of  tne  data  is  large.  Menael"'’  g-^ves  tne  ''ule  o*  thumt 
that  tne  coe-^ficient  of  variation  should  not  exceed  lOL. 

In  linear  regression  J  data  values  are  least  squares  fitted  to  a  straight 
line  wncse  equation  is  given  by; 

X.I  -  a.  7k  -t-  I  /. 


where  y  is  the  dependent  variable,  x  is  the  indeoendent  variable,  a  is  the  slope 
and  b  is  the  intercept.  Assuming  that  all  the  inaccuracies  occur  in  the 
measured  values  of  y  and  that  tne  parent  standarc  deviation  of  the  data,  r, 
is  independent  of  x,  the  variances  and  covariance  of  the  parameters  a  and  b 
can  be  expressed  as : ' ^ 

\  /  r,  ■ 

V  :  '''■  j 


I'X; 


V(bi  =  S-' 

A  lx) 


-CV  t. 


'T  0; 


where  i  (x)  =  J 
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46. 


j, 

K^  f/u-xV^  _  i,/.  ^ 


(A53) 


(A54) 


In  Eqs.  (A52)-(A54)  it  is  assumed  that  the  covariances  between  a^,  a^  and 

A 

In  S  are  zero  and  that  the  number  of  inert  strength  samples  is  the  same  as  the 
number  of  samples  tested  in  fatigue  at  each  applied  stress. 
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